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Abstract: We present the electroweak spectrum for the Next-to-Minimal Super- 



> 

O symmetric Standard Model at the one-loop level, e.g. the masses of Higgs bosons, 

sleptons, charginos and neutralinos. For the numerical evaluation we present a 

(^ mSUGRA variant with non-universal Higgs mass parameters squared and we com- 

^ pare our results with existing ones in the literature. Moreover, we briefly discuss the 

r| implications of our results for the calculation of the relic density. 
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1. Introduction 

Supersymmetric extensions of the standard model (SM) are promising candidates for 
new physics at the TeV scale [1, 2, 3] as the solve several short-comings of the Stan- 
dard Model (SM). The Minimal Supersymmetric Standard Model (MSSM) solves 
the hierarchy problem of the SM [4], leads to a unification of the gauge couplings [5] 
and introduces several candidates for dark matter depending on how SUSY is broken 
[6, 7]. On the other hand, a new problem arises in the MSSM: the superpotential 
contains a parameter with mass dimension, namely the so called /i parameter which 
gives mass to the Higgs bosons and higgsinos. From a purely theoretical point of view, 
the value of this parameter is expected to be either of the order of the GUT/Planck 
scale or exactly zero, if it is protected by a symmetry. For phenomenological aspects, 
however, it is necessary that it is of the order of the scale of electroweak symmetry 
breaking (EWSB) and it has to be non-zero to be consistent with experimental data. 
This discrepancy is the so called yU-problem of the MSSM [8]. 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [9] provides 
an elegant solution to this problem. The particle content of the MSSM is extended 
by an additional gauge singlet, which receives a vacuum expectation value when 
supersymmetry is broken. The corresponding term in the superpotential gives then 
rise to an effective /x-term which is naturally of the order of the EWSB scale. Also 
in this model several regions exist in parameter space where one obtains the correct 
relic density to explained the observed dark matter [10, 11]. It turns out that in 
high scale models like mSUGRA several regions exist which are rather sensitive to 
mass differences of the various supersymmetric particles, in particular the masses 
of the Higgs bosons, the neutralinos and the staus, the supersymmetric partners of 
the tau-lepton, and require precise calculations of these masses. The corresponding 
regions are the so-called Higgs funnel(s) and the co-annihilation regions. Motivated 
by this observation we calculate the Higgs masses, the neutralino masses and the 
staus at the one-loop level. 

This paper is organized as follows. We first detail our calculation of the mass 
spectrum in sec. 2. In sec. 3 we present the constrained NMSSM, which serves as our 
reference scenario and perform a numerical analysis of our implementation, sec. 4 is 
devoted to a comparison of our results with the public program package NMSSM-Tools 
[12]. Finally, we give a few examples for the calculation of the dark matter relic 
density in sec. 5 and draw our conclusions in sec. 6. We collect the couplings and 
one-loop self-energies in the appendix where we include for completeness also those 
for the Z-boson and neutral Higgs bosons which have already been given in ref . [13] . 

2. Calculation of the One-Loop Mass Spectrum 

In this section we fix our notation and discuss briefiy the DR renormalization of the 



relevant masses, where we follow closely ref. [14]. 

2.1 Superpotential and soft SUSY breaking terms of the NMSSM 

As already stated above, the solution to the /i-problem of the MSSM is the replace- 
ment of the bilinear /x-term by a coupling between the Higgs superfields and an 
additional gauge singlet S leading to the superpotential 

W^NMSSM = -HuQYuU' + HdQYdb' + HdLY.E'' + XHuH^S + ^-kSSS. (2.1) 

where the last term is introduced to forbid a Peccei-Quinn symmetry which would 
lead to an axion in contradiction to experimental results, see e.g. ref. [15] and refs. 
therein. Moreover, we have only taken into account dimensionless couplings to avoid 
the /i-problem of the MSSM. 

The scalar component S oi S receives after SUSY breaking a vacuum expectation 
value (VEV), denoted Vg, which leads to 

/Xcff = -J^^'"s, (2.2) 

where we have used the decomposition 

5 = -^ (0. + leys + Vs) . (2.3) 

Since Vg and thus also /icg are a consequence of SUSY breaking one finds that /Xes is 
naturally of the order of the SUSY breaking scale. 

All interactions are fixed by the gauge structure and the above superpotential. 
We have used the Mathemtica package SARAH [16] to calculate all vertices, mass 
matrices including the one-loop corrections and renormalization group equations of 
the model. 

In the following, we use the standard conventions, where for a matter superfield 
X, X denotes its scalar component and X denotes its fermionic component. In case 
of the Higgs fields and the gauge singlet, Hu^/S are the scalar components, while 
Hu,d/S are the fermionic higgsinos and the singlino. 

At tree level the scalar potential receives contributions from several sources: from 
the the superpotential in eq. (2.1) the so-called F-terms given by 



vp = E 



dwi<p,) 



d(f)i 



(2.4) 



The sum runs over all chiral superfields (pi, which are then replaced by their scalar 
component (pj. The D-terms are 

^D = ^EE|E^*^.>^-r' (2-5) 

a a i,j 



and finally the soft breaking terms 

VsB,2 = m^„ \Hu? + mlji/rfp + m||5|2 + gt^|g + 
+Vm\L + !)^m\b + U^mlil + 

+ ^ (Mi BB + M2 iy„iy" + M3 ^,^" + /i.e.) (2.6) 

VsB,3 = -H^QTuU^ + H^QTdD^ + Ham& + TxHuH^S + ^T^SSS (2.7) 

The sum in eq. (2.5) runs over all gauge groups g and over the corresponding gen- 
erators a, i.e. ^A" in the case of SU{3), ^a" in the case of SU(2), and |F^ for the 
U{1). Here, A" are the Gell-Mann matrices, a" the Pauli-matrices, and Y is the 
hyper charge. 

2.2 Minimum Conditions of the Vacuum 

Once electroweak symmetry gets broken, both Higgs doublets receive a VEV and we 
decompose the scalars similar to eq. (2.3) 

Hu,d = -y= {4>u,d + icru4 + Vu,d) ■ (2.8) 



At tree level, the minimum conditions for the vacuum are the so-called tadpole 
equations 

(2.9) 



Ti = 

dvi 



=0,o-i=0 



with 
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--vlvu^e{K\] —VsVu^e{Tx], (2.10) 

-^r;2T;^Re{ft:A} - -^VsVdRe{Tx], (2.11) 

Tg = — = m|t;s + t;^|/€|^ - VdVsVu^e{K.>^] + 2 (^rf + ^«)^s|Ap 

+ -^t;2Re{T4 - l^VdVu^e{Tx]. (2.12) 

Here we have chosen a phase convention where all VEVs are real. For the later 
calculation of the one-loop corrections to the Higgs boson masses one needs the 
evaluation of the tadpole equations at the one- loop level, leading to corrections 5ti . 
As renormalization condition we demand 

Ti + 5tf^ = for i = rf, M, s. (2.13) 



All calculations are performed in 't Hooft gauge using the diagrammatic approach. 
The explicit formulas for 5t\ are given in app. D. In our subsequent analysis we will 
solve eqs. (2.13) for the soft SUSY breaking masses squared: m'j^^,m'jj^, and 
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All parameters in eqs. (2.10)-(2.12) are understood as running parameters at a 
given renormalization scale Q. Note that the VEVs Vd and Vu are obtained from the 
running mass mz{Q) of the Z-boson, which is related to the pole mass mz through 



m|(Q) 



91 + 92 
4 



{vl + vl)=ml + Re{Yilz{ml)]. 



(2.14) 



The transverse self-energy 11^^ is given in app. E.l. Details on the calculation of 
the running gauge couplings at Q = mz can be found in ref. [14]. The ratio of these 
VEVs is denoted as in the MSSM by tan/3 = Vu/vd- 



2.3 Masses of the Higgs bosons 

The tree-level mass matrices for the neutral scalar Higgs bosons and pseudo scalar 
Higgs bosons can be calculated from the scalar potential according to 
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(2.15) 



Afe=0,o-fc=0 



respectively, with i,j = 1,2,3 = u, d, s. The matrices are symmetric and the entries 
in case of the scalar Higgs bosons are 
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T,33 
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+ St^.Vr + ^(^^ + vl) |A|2 + V2vsRe{T^] - VdVuRe{K\}, (2.21) 



while those of the pseudo-scalar ones are 

^T,n =<, + \ [vl + ^s) |AP + ^ {gl + gl) {vl - v^:) , (2.22) 



Id 2 
m;;,'Vo = -^v,Re{Tx] + -t;,^Re{/tA}, (2.23) 



4'i° = -7fs^(i{Tx] + -r^i 
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ml. + M + ^s) |AP + « (^? + gl) {vl - ^d), (2.25) 



2 



"^T23 = ^^^^^^JTa} -t;dt;sRe{KA}, (2.26) 

v2 

^T,33 = ^1 + ^''^^ + ^ (^^ + ^«) 1^1' - v^^^ReJTj + t;,t;„Re{ft:A}, (2.27) 

where rnl^, m|^^ and m| satisfy the tadpole equations. 

The diagonalization of the mass matrices rrir^ and irir^ leads in total to five 
physical mass eigenstates and one neutral Goldstone boson which becomes the lon- 
gitudinal component of the Z-boson. The five physical degrees of freedom are: three 
CP-even Higgs bosons denoted /ii,2,3 and two CP-odd bosons denoted A? 2- The 
corresponding rotation matrices Z'^ und Z"^ are defined through 

^x^2,x^x,T ^ ^2,x^^ X = /l, A° . (2.28) 

Moreover, we note that we order all masses such, that nii < rrij Hi < j. 

The one-loop scalar Higgs masses are then calculated by taking the real part of 
the poles of the corresponding propagator matrices 



Det 

where 



2,h, 



ptl-mttip') =0, (2.29) 



mli{p')=m'/-U,,{p'). (2.30) 

Here, rhh is the tree-level mass matrix from eq. (2.15). Equation (2.29) has to be 
solved for each eigenvalue p^ = mf. The same procedure is also applied for the 
pseudo scalar Higgs bosons. The complete 1-loop expressions for the self energy of 
the CP-odd and even Higgs bosons are given in apps. E.2 and E.3. 

The charged Higgs sector consists of H^ and H:^. The mass matrix in the basis 
(^H^ , H^'*) is diagonalized by a unitary matrix Z^ 

ZW'^^Z+'t = ^^;f;. (2.31) 

The eigenstates yield as in the MSSM the longitudinal component of the W^-boson 
and a charged Higgs boson H~^ with mass 

^2,H^ ^ jvl + vl)i2vsRe{K\} + v^v^jgl - 2|Ap) + 2V2v,Re{Tx}) 



and the one-loop mass 



m 



IL 



rrir, 



2,H+ 



-Re{UH+H+{ml'L )} 



(2.33) 



where the self-energy can be found in app. E.4 

2.4 Chargino and neutralino masses 

As for the Higgs bosons discussed in the previous section, one has to find the real 
parts of the poles of the propagator matrix to obtain the masses of charginos and 
neutralinos. At the tree-level the chargino mass matrix in the basis ip" = {W^ , H^)^, 
ip+ = (iy+, iJ+) is given by 



£ r+ = -i!-^Mf^+ + h.c. 



(2.34) 



with 



Mf 




(2.35) 



This mass matrix is diagonalized by a biunitary transformation such that U*M^ V'^ 
is diagonal. The matrices f/ and V are obtained by diagonalizing M^ (M^ )''' and 
(M^ )*(My )^, respectively. At the one-loop level, one has to add the self-energies 



Ml^ (pI 



Mf 



^tip')-^Upt)M} -M^ ^Uvi). 



(2.36) 



In case of the neutralinos one has a complex symmetric 5x5 mass matrix which 
in the basis ip^ = {B, W^, H^, if °, S)'^ is at the tree-level given by 
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\g\Vu -\92Vu 





-\9\'Vd \g\Vu 
\92vd -\92VU 

ynVuX -^VdX V2vsK J 



(2.37) 
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One can show that for real parameters the matrix M^ can be diagonalized 
directly by a 5 x 5 mixing matrix A^ such that N*M^ N"^ is diagonal. In the complex 
case, one has to diagonalize M^ (M^ )'''. At the one-loop level we obtain 



mUpI 



M^-^ 



+Mf (s°/(pD + S° (K 



Mf 



(2.38) 



The complete self-energies for neutralinos and charginos are given in apps. E.5 and 
E.6, respectively. 



2.5 Masses of sleptons 

In the basis {ei, fiL,TL,eji,fiji,fR), the mass matrix of the charged sleptons at the 
tree-level is given by 

with the diagonal entries 

mlL = ml + |(Fe)*(V;)^ + l{gl- gj) [vj - vl)l„ (2.40) 

mln = -I + %Y,f{Y^r + f (^^^ - ^^^jls- (2.41) 

Where I3 is the 3x3 unit matrix. This matrix is diagonalized by a unitary mixing 
matrix Z^: 

ZV^^Z^^ = rr^^ . (2.42) 

The corresponding mass matrix at the one-loop level is again obtained by taking into 
account the self-energy according to 

and the one-loop masses are obtained by calculating the real parts of the poles of 
the propagator matrix. The expression for Ilij{pf) can be found in app. E.7. 

Finally, in the basis [i)e,i'fj,,i'r) the tree- level sneutrino mass matrix is given by 

mr = ml + l[gl + gl)[vl-vl)u (2.44) 

This matrix is diagonalized by a unitary mixing matrix Z^: 

Z^m'/Z^^ = rr^li^ . (2.45) 

Similarly as above the one loop mass matrix is given by 

ml'lip^)=m'/-U,,ip^). (2.46) 

The one-loop masses are obtained by calculating the real parts of the poles of the 
propagator matrix. The expression for n,>;>(pf) can be found in app. E.8. 

3. The constrained NMSSM 

3.1 The model and its free parameters 

In the subsequent numerical analysis, we are mainly interested in precision calculation 
of the SUSY masses and potential effects in the calculation of the relic density. To 



reduce the number of free parameters we therefore focus on a scenario motivated 
by minimal supergravity (mSUGRA) [17]. More precisely, we study a variant of the 
constrained NMSSM [18, 19] where we allow for non-universal Higgs mass parameters 
squared at the GUT scale. In our setup, these parameters are determined with the 
help of the tadpole equations (2.13) at the electroweak scale. As a side remark, 
we note that also other recently used mSUGRA versions of the NMSSM contained 
non-minimal features either for the scalar mass parameter and/or for the trilinear 
couplings. 

We apply the following boundary conditions for the gaugino masses Mi, M2, M3 
and the soft breaking masses of the squarks and sleptons mf at the GUT scale, which 

is defined as the scale where the f/y(l) and SU{2)i couplings fulfill \ \gi = 92- 



Ml = M2 = M3 = Ml/2, (3.1) 

2 2 2 2 2 2-1 /oo\ 

m^ = rriQ = m^ = m^ = m^ = tUq I3. (3.2) 

The trilinear scalar couplings Tj are given by 

T„ = AqYu, Td = AoYd, Te = AoYe, Tx = AxX, and T^ = A^k. (3.3) 

Here, Aq is defined at the GUT scale, while A, k, Ax and A,^ can be defined either 
at the GUT or at the SUSY scale. Together with the values for tan/3 = — and f, 
the spectrum is fixed. To summarize, we have nine input parameters, 

Ml/2, fuo, Aq, a, k. Ax, A^, Vs, and tan/3. (3.4) 

Note, that we allow for non-universalities in the trilinear parameters for an eas- 
ier comparison with the existing literature but in principal we could take all A- 
parameters equal at the GUT scale. We choose in the following f ^ > and A, /t G 

[-1,1]- 

3.2 Procedure to evaluate the SUSY parameters at the electroweak scale 

In order to connect the parameters at various scales, we use the renormalization 
group equations (RGEs), which are calculated at the two-loop level in the most 
general form using the generic formulas given in ref. [20]. We have compared the 
obtained expressions for the RGEs with those given in ref. [15] in the limit where 
only the third generation Yukawa couplings contributes. There has been a slight 
difference in the two-loop /3-function of Ax = Tx/X, but it was confirmed by the 
authors of ref. [15] that our result is correct. The RGEs themselves can easily be 
calculated by the CalcRGEs command of SARAH and a print-out can be found at [21]. 
In the calculation of the gauge and Yukawa couplings we follow closely the proce- 
dure described in ref. [22] : the values for the Yukawa couplings giving mass to the SM 
fermions and the gauge couplings are determined at the scale Mz based on the mea- 
sured values for the quark, lepton and vector boson masses as well as for the gauge 



couplings. Here, we have included the one-loop corrections to the mass of W- and Z- 
boson as well as the SUSY contributions to 5vb for calculating the gauge couplings. 
Similarly, we have included the complete one-loop corrections to the self-energies 
of SM fermions extending the formulas of [14] to include the additional neutralino 
and Higgs bosons. Moreover, we have resummed the tan/3 enhanced terms for the 
calculation of the Yukawa couplings of the 6-quark and the r-lepton as in [22] . The 
vacuum expectation values Vd and Vu are calculated with respect to the given value 
of tan/3 at Mz- Furthermore, we solve the tadpole equations to get initial values 
for m?^ , rn^j and m|. Afterwards the RGEs are used to obtain the values at the 
GUT scale and all boundary conditions including A and k are set as described above. 
Then, an ROE running to the SUSY scale is performed and the SUSY masses are 
calculated at the one-loop level and for the neutral and pseudo scalar Higgs bosons 
we include beside the one-loop contributions presented here also the known two-loop 
contributions [13]. For this purpose also the numerical the values for the VEVs at 
MsusY are needed. These are derived using the two-loop RGEs 

P.. = -v^ [if^ + iT) (3.5) 

with i = u,d. Here, 7} and 7} are the anomalous dimensions for the two Higgs- 
doublets at the one- and two-loop level, respectively. The corresponding expressions 
are given in app. B. Let us recall that the input value for Vg is already given at 
^susY- These steps are iterated until the masses converge with a relative precision 
of 10^^. The complete procedure has been implemented in SPheno [22]^. 

3.3 An example spectrum 

In Table 1 we give as an example the masses of the Higgs bosons, chargino, neutralinos 
and third generation sfermions at tree-level as well as at the one- and two- loop level 
for the parameter set 

mo = 180 GeV,mi/2 = 500 GeV, Ao = A^ut = _i5oo GeV, A^^^ = -36 GeV, 
tan/3 = 10, K^UT _ g.ll, A^^t _ g.l, vs = 13689 GeV . 

which is close to the benchmark scenario 1 of ref. [18]. As can be seen in Table 
1, the corrections are sizable ranging from 0.1 % to 23.6 % in case of the lightest 
Higgs boson. This large correction is well known and the main reason for including 
the two-loop corrections. The corresponding two-loop Higgs masses as well as the 
relative correction with respect to the one- loop results are also displayed in Table 1. 
Again the largest correction with 5.2 % is in case of the lightest Higgs boson mass. 
As an estimate of the remaining theoretical uncertainty we have varied the renor- 
malization scale in SPheno. We show in fig. 1 the scale dependence for masses of 
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Particle 


iriT [GeV] 


miL [GeV] 


A [%] 


m2L [GeV 


A [%] 


hi 


86.7 


113.3 


23.5 


119.6 


5.2 


h2 


863.1 


934.2 


7.6 


937.3 


0.3 


hs 


2073.9 


2073.9 


< 0.1 


2073.9 


< 0.1 


A? 


76.4 


69.3 


10.2 


69.5 


0.3 


AO 


865.2 


937.2 


7.7 


940.4 


0.3 


X? 


211.6 


207.6 


1.9 


- 


- 


X^ 


388.2 


398.4 


2.6 


- 


- 


xl 


987.9 


980.5 


0.7 


- 


- 


xl 


993.0 


985.1 


0.8 


- 


- 


xl 


2074.8 


2074.9 


< 0.1 


- 


- 


xt 


388.2 


398.6 


2.6 


- 


- 


xt 


993.3 


985.9 


0.7 


- 


- 


Tl 


191.1 


193.3 


1.2 


- 


- 


T2 


388.1 


393.1 


1.1 


- 


- 


tl 


506.9 


541.8 


6.4 


- 


- 


k 


914.4 


949.3 


3.7 


- 


- 


hi 


845.3 


880.4 


3.9 


- 


- 


&2 


961.9 


1008.5 


4.6 


- 


- 


^ 


1107.2 


1154.2 


4.1 


- 


- 



Table 1: Comparison of the tree- level rriT and loop masses at 1-loop {miL) and 2-loop 
(?712l). a is the relative difference |1 — ^^^^| respectively |1 



miL I 



neutral scalar Higgs boson at the one- and two-loop masses normalized to their values 
at Q = 1 TeV and vary the renormalization scale Q between 200 GeV and 2.2 TeV. 
As can be seen, the large variation of 8% at one-loop for the lightest Higgs, which 
is mainly the lighter SU{2) doublet Higgs in this case, is reduced at two-loop to less 
than 2%. In case of the heavier Higgs bosons the scale dependence is significantly 
smaller showing a significant improvement when going from the one-loop level to the 
two-loop level. However, we remark that the values of A and k are small in this 
scenario and we expect a stronger dependence in case of larger couplings. 

The picture changes slightly in the case of the pseudo scalar bosons as can be 
seen in fig. 2. While the heavier pseudo scalar behaves exactly as the second scalar 
field since both originate to 99.5% from if^, the scale dependence for the lighter 
pseudo scalar is smaller compared to the lightest scalar field, but hardly improves 
at the two-loop level. This is because in the two-loop part contain 'only' the strong 
contributions of the third generation squarks whereas this state is mainly a singlet 
state and, thus, the contributions due to the Yukawa couplings would be needed for 
a further improvement. 
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Figure 1: Dependence of CP even Higgs masses on the renormalization scale Q at 1-loop 
(red) and 2-loop level (dashed blue) normalized to the value at Q = 1 TeV. From left to 
right: m^.^, ruh^ and mh^. 
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Figure 2: Dependence of CP odd Higgs masses on the renormalization scale Q at 1-loop 
(red) and 2-loop level (dashed blue) normalized to the value at Q = 1 TeV. Left: m^o. 
Right: TTT-^o. 

In Figure 3 the scale dependence for different neutralinos is shown. As can be 
seen, in case of the three lighter states the scale dependence is reduced from the level 
of about 1.5% to 3-5 per-mill. In case of the singlet state X5 the scale dependence is 
already small due to the small values of A and k. We note that the scale dependence 
of Xi ixt ^^^ Xi) is nearly the same as that of X2 (Xs) as these state have their 
main origin in the same electroweak multiplet. 

Finally we show in fig. 4 the scale dependence of the staus. The scale dependence 
at tree level amounts to about 2-2.5% and is reduced at one-loop level to about 1% 
and less where the fi shows still the larger dependence. The sleptons of the first two 
generations show a somewhat smaller scale dependence as in their cases the Yukawa 
couplings do not play any role. 

4. Comparison with the hterature 

To date, the program package NMSSM-Tools [12] has been the only complete spec- 
trum calculator for the NMSSM. NMSSM-Tools uses for the constrained NMSSM the 
parameters mo, M1/2, ^0 and A^ at the GUT scale whereas tan/3 and A are given at 
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Figure 3: Dependence of the masses of the hght neutrahnos on the renormahzation scale 
Q at tree (red) and 1-loop level (dashed blue) normalized to the value at Q = 1 TeV: From 
left to right and from above to below: m,-o, m^o, m^a and m^a. 
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Figure 4: Dependence of the stau masses of the sleptons on the renormahzation scale Q 
at tree (red) and 1-loop level (dashed blue) normalized to the value at Q = 1 TeV. 

the electroweak scale. Moreover, in NMSSM-Tools the tadpole equations are solved 
with respect to \vs\, k, and m|. We have performed a detailed numerical comparison 
of our implementation with the version 2.3.1 and present here a few typical examples. 

4.1 Differences between the programs 

Since both programs use different methods to calculate the spectrum, we have also 
done a comparison where we modified the codes such that both codes use equivalent 
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methods except for small details. First, the implementation of NMSSM-Tools involves 
two different scales, namely the SUSY scale defined as 

QlvSY = MksY = \ (2<, + <n + ^L) ' (4.1) 

and the scale at which the masses are calculated, 

QsTSB = f^qif^m- (4.2) 

In SPheno, all masses are evaluated at the SUSY scale, so that we had to set 
QsTSB = QsusY in the relevant routines of NMSSM-Tools. Second, as already stated 
in sec. 3.2, the two-loop /3 function of A\ has been corrected in the public version 
of NMSSM-Tools. However, in general the numerical effect on the spectrum is rather 
small. 

In the Higgs sector the loop contributions are taken into account differently 
in both codes. While SPheno takes the complete one-loop correction including the 
dependence of the external momenta, NMSSM-Tools uses the effective potential ap- 
proach, e.g. setting the external momenta to zero. NMSSM-Tools calculates afterwards 
the momentum dependent contributions from top and bottom quarks. Also the in- 
cluded contributions differ: in SPheno the complete one-loop corrections to both, 
scalar and pseudo scalar Higgs bosons, and the two-loop contributions as given in 
ref. [13] are included. In NMSSM-Tools beside the dominant contributions due to 
third generation sfermions also electroweak corrections and some leading two-loop 
corrections for the scalars are calculated, while for the pseudo-scalars only the dom- 
inant one-loop corrections due to tops, stops, bottoms, and sbottoms are included. 
In addition, some corrections due to charginos and neutralinos are absorbed in a 
redefined Ax. To account for these differences we have switched off the two- loop 
parts in both codes. Furthermore, we have set the external momenta of the loop- 
diagrams of scalars in SPheno to zero. Finally, we have kept only those corrections 
to the pseudo-scalar masses in SPheno which are also included in NMSSM-Tools, but 
neglected the additional corrections absorbed in Ax- In the following, we refer to 
these modified versions by SPheno mod and NMSSM-Tools mod, respectively. 

Also in the chargino and neutralino sector the implementations are different: in 
SPheno the complete one-loop corrections are implemented whereas in NMSSM-Tools 
the corrections to the parameters Mi, M2, and /iefi are taken into account. In the 
slepton sector the differences are largest: SPheno contains the complete one- loop 
corrections whereas in NMSSM-Tools the calculation is done at tree-level. Last but 
not least we note that the data transfer has been done using the SLHA2 conventions 
[23]. 

4.2 Results of the comparison 

As a first reference scenario, we take the benchmark point 1 proposed in ref. [18]. 
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The corresponding input parameters for NMSSM-Tools are 

mo = 180 GeV, mi/2 = 500 GeV, Aq = -1500 GeV, tan/3 = 10, 

A^usY = 0.1, A^UT = _33.45^ ^^^ > q. (4.3) 

In the following we will vary rriQ and the keep the other parameters to the values 
shown here. 

In the left graph of fig. 5, we show the mass of the lightest scalar hi as a function 
of rriQ. The largest discrepancies arise for the lighter scalar and pseudo scalar boson, 
where the relative differences between the complete calculation of both programs 
amount up to 2.5 and 35%, respectively. In case of h^ this is a combination of 
the p^ terms in the loop-functions and the additional two-loop contributions. The 
differences in case of A^ can easily be understood by noting that in NMSSM-Tools 
only the contribution of third-generation sfermions are taken into account whereas 
we include the complete one-loop corrections plus the known two-loop contributions. 
In case of the modified program codes these differences reduce to at most 2% which is 
meanly due to two differences: (i) the way the top Yukawa coupling is calculated and 
(ii) the way the tadpole equations are solved which leads to somewhat different values 
between the two programs. There is no visible difference between NMSSM-Tools and 
NMSSM-Tools mod for the pseudo scalar and the heavy scalars. The reason is that in 
the case of the pseudo scalar no two-loop corrections are calculated in NMSSM-Tools 
and in case of the heavy scalars they are very small. 

Finally, we have also cross-checked our results in the Higgs sector with ref. [13] 
and we have found agreement better than one per-mill when using the set of soft 
SUSY parameters at the scale Qstsb- This small difference is an effect of the Yukawa 
and scalar-trilinear couplings of the first two generations which we take also into ac- 
count. If we restrict ourself to third generation couplings there is an exact agreement 
between both calculations. 

Concerning the chargino and neutralino masses, the agreement between the two 
spectrum calculators is rather good as can be seen in fig. 6. The relative differences 
are at most 1% and in general slightly below 0.5%. In case of the sleptons the 
differences are more pronounced as can be seen in fig. 7 which is due to the differences 
between tree-level and one-loop calculation and amounts in up to 3% and 0.6% for the 
light and heavy stau, respectively. Note, that although for LHC physics one expects 
similar experimental uncertainties, this precision necessary for a future linear collider 
or dark matter calculation require the inclusion of the radiative corrections to the 
slepton masses. 

5. Effects on the relic density of dark matter 

It is well known that the prediction of the dark matter relic density ^cdm^^ is very 
sensitive to the exact mass configuration of the scenario under consideration [24] . For 
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Figure 5: Comparison of the masses in GeV of the hghtest scalar (upper left) , the 
lightest pseudo scalar (upper right), heavier scalar masses (lower left) and heavier pseudo 
scalar mass (lower right) as a function of tjiq (in GeV). All other parameters are fixed 
as in eq. (4.3). The lines correspond to: are for unmodified version of SPheno (full red), 
NMSSM-Tools (dotted blue), SPheno mod (dashed red) and NMSSM-Tools mod (dot-dashed 
blue). 



a neutralino LSP, this is, e.g., the case for the annihilation through Higgs- resonances, 
but also in the case of neutralino-sfermion co-annihilation. For the latter, the mass 
difference between the two particles plays a key role in the calculation of the resulting 
relic density. Therefore, it is necessary to calculate the complete spectrum as precisely 
as possible to get viable results of allowed regions of parameter space with respect 
to the constraints imposed by the presence of dark matter. Let us recall that recent 
measurements by the WMAP satellite in combination with further cosmological data 
lead to the favored interval 



0.1018 < ncDufi^ < 0.1228 



(5.1) 



at 3o" confidence level [25]. 

We compute the relic density of the lightest neutralino using the public program 
package mi crOMEGAs 2.4.0 [26]. To this end, we have implemented the NMSSM 
particle content and corresponding interactions into a model file for CalcHEP [27], 
which is used by micrOMEGAs to evaluate the (co) annihilation cross-section. The 
relevant interactions have again been calculated and written into the model files by 
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Figure 6: Comparison of cliargino and neutralino masses (in GeV) as a function of tuq 
(in GeV). All other parameters are fixed as in eq. (4.3). The lines correspond to the 
unmodified versions of SPheno (full red) and NMSSM-Tools (dashed blue). Up left: light 
neutralinos Xi- Up right: neutralino X2 and chargino xf (SPheno: black dotdashed, 
NMSSM-Tools: black dotted). Down left: neutralinos X3 (thin lines), X4 (thick lines) and 
chargino xt (SPheno: black dotdashed, NMSSM-Tools: green dotted). Down right: X5- 

the program package SARAH. Let us note, that we take into account important QCD 
effects, such as the running strong couphng constant and the running quark masses 
[28, 29, 30]. 

As an example, we illustrate the effect of the one-loop correction to the slepton 
masses on the dark matter relic density in a region of dominant neutralino-stau 
coannihilations. In fig. 8, we show the isolines corresponding to the upper and lower 
limit of eq. (5.1) in the mQ-mi/2 plane. All remaining parameters of eq. (3.4) are 
fixed to 



15, K 



^GUT 



tan/3 
= 30 GeV, Ao 



SUSY 



^GUT 



-0.05, ASUSY = _o.l, 

= 1000 GeV, Vs = 2- 10^ GeV. (5.2) 



One clearly sees that the allowed parameter range gets shifted depending on the 
precision with which the spectrum is calculated. More, the two regions shown do not 
overlap as can also be clearly be seen in the figure to the right. 
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Figure 7: Comparison of selectron and stau masses (in GeV) as a function of mo. All 
other parameters are fixed as in eq. (4.3). The lines correspond to the unmodified versions 
of SPheno (full red) and NMSSM-Tools (dashed blue). 
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Figure 8: The isolines corresponding to JIcdm^^ = 0.1018 and QcDuh'^ = 0.1228 in 
the mo~mi/2 plane for dominant neutralino-stau coannihilations. All other parameters 
are fixed as in eq. (5.2). The red solid lines have been obtained for the complete mass 
spectrum at the one-loop level, while for the black dashed line the loop corrections to the 
slepton masses have been disabled. The right graph corresponds to a zoom into the left 
one. 



For a point with ^cdm^^ = 0.112 at mi/2 — 451.2 GeV, the resuhing one-loop 
corrected masses of the hghtest neutrahno and the hghter stau are tti^o = 186.0 GeV 
and r/ifj = 196.8 GeV, respectively. In consequence, coannihilations account for 
about 60% of the total annihilation cross-section, where the most important final 
states are rhi (27%) and tZ^ (15%). A sizable contribution of about 14% (5%) 
comes also from stau-antistau (stau-stau) annihilation. The remaining contributions 
are mainly from neutralino pair annihilation. For lower values of mi/2 ^ 200 GeV, 
the coannihilations become less important within the WMAP-favored region, the 
dominant mechanism is then neutralino pair annihilation into r^r^ pairs through 
stau-exchange. 
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6. Conclusion 

The NMSSM is an attractive extension of the MSSM, in particular as it solves the 
/i-problem of the MSSM and as it leads to new phenomenology at present and future 
collider experiments. It can also explain the observed amount of dark matter in 
the universe. However, in particular for comparison of the WMAP data improved 
theoretical predictions are necessary. We therefore have presented in this paper the 
complete one-loop calculation of the electroweak sector: Higgs bosons, charginos, 
neutralinos and sleptons. While in case of the Higgs bosons we have reproduced 
known results, the corrections to the other particles have not yet been discussed in 
the literature. We have shown that the corrections amount to the order of a few per- 
cent. While the corrections are most likely below the precision of the coming LHC 
data, they are clearly important for comparison with WMAP data and also with a 
future international linear collider, and thus crucial for precision investigations of the 
NMSSM parameter space. 
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A. Squark mass matrices 

For completeness, we display here the mass matrices of the squarks, since they also 
enter the one-loop corrections and are particularly important in the case of the 
Higgs-bosons. In the basis (c^l, sl, 6l, (Ir, sr, BrJ , the mass matrix for the down-type 
squarks is given by 



m. 




2,d 

2,d I m 

T 



- VsVu\Y2) rnA ] ' 



Iv^Tl - VsV^XY*) 



where the diagonal entries read 






•'LL — "•'Q ^ 2 '^ '^ 24 

12 



^ + gl 1^,2 

2 



^'A = < + ^yJy: + '-^M- -i) 13- (a.2) 
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The corresponding expressions for the up-type squarks in the basis (m^, c^, t/,, ur, cr, t^) 



are 



with 






^KR 



2 I ^w-iz*-izT I '^fi'2 9l / 2 2\-i 



2 2 

,2,1* _ ^2 I "^n-i^T-i^* I 9l / ,2 „,2 



These matrix are diagonahzed by unitary mixing matrices Z'^: 

Z^m^^Z^^ = m^f g , q = d,u. (A.5) 

B. Anomalous Dimensions 

In this app., we give the detailed expressions of the anomalous dimension of the 
Higgs-fields, which are needed for the RGE evaluation of the VEVs. 

7^^^ = 3Tr(F,Fi) - ^gf - \gl + |Ap + Tr(y;i;t) (B.l) 

+ 16(?3'Tr(r,ri) + yiTv{Y,Y}) - 3|A|2Tr(KFj) - 9Tr(F,Fjr,r]) 
- 3Tr(F,Fjy;F„t) - 3Tr(F,y;ty;v;) (B.2) 

I'l = 3Tr(y;y;t) - ^^gl - Igl + \X\' (B.3) 



n: = ^91 + ^9W2 + ^^2 - 2|A| VP - 3|Ar - 3|ApTr(r,r,' 



|A|^Tr(y;r;) + -g^TrlY^Y^^ + IGgs^TrlYX) ' STrlYdYjYX] 
9Trfy;y;ty^yA (B.4) 



7f) = 2|/s:|2 + 2|Ap (B.5) 

7f = l9l\M' + 6^2'lAr - 8|«:r - 8|A| Vr - 4|Ar - 6|ArTr(r,ri) 

- 2|A|2Tr(FeF;) - 6\\\^Tt(yX) (B.6) 

C. Couplings 

We list in the following all couplings of the NMSSM which contribute to the elec- 
troweak self-energies or influence the annihilation or coannihilation of the neutralino. 



19 



These and all other couplings of the NMSSM can be derived with the command 
MakeVertexList [EWSB] of SARAH. A pdf version is also available at [21]. 
We define the following abbreviations: 

\^+gl^gl-A\^ (C.l) 

X^gl- 2\' (C.2) 

gl = 9l - 9l (C.3) 

9l^9l + 9l (C.4) 

Ai = 2v,k\* + 2vs\iC + ^2 2 Re{TA} (C.5) 

As = -2vu\ + VdK (C.6) 

K^ = -2vd\ + v^K (C.7) 

Furthermore, ce is cos{Qw) and sq is sin(0i4/). 

C.l Two fermion - One Scalar 

^t?,xi>.3 =l{- 92NI2KA + v^AAT* sAT^.Zj, + V2\Nl,Nl,Z^^, 

- 9iNUNl,ZZ, + g^Nl.Nl.Zl^^, + v^AiV^^iV* 3ZJ2 + 92Nl,Nl,Z^^^ 
+ 9iNl, {KsZI^^, - Nl,ZZ,) + V2\Nl,Nl,ZZ, - 2V2kNI,NI,Z^^, 

+ Nl,{g,Nl,Z^^, - g^Nl.ZZ, + v^A(Ar*4Z,^3 + iV^s^^))) 

(C.8) 

r|iX?,/.*3 = I {^"' (^*^' (^'^*^' " ^'^*^0 + ^1^*11^*^3 - ^2Ari,2Art,3 + V2X*Nt,,Nt 
+ v^A* Ari^4A^t,5) + v^^tfs ( - S/t^AT.^sATi^s + A* (^Nt.^Nt,, + N^.^Nt, 
+ ZJ2 (iVi,4 ( - (7iiVt2i + 92Nt,2) + ( - 9iNt,i + 92Nt,2)Nt,4 
+ V2\* (Ar,^3iVt25 + Nt,,Nt,s) ) ) (C.9) 

- ^lA^^iV^i^i + 92K,K,Z^^, - V2\Nl,Nl,Zf^, + 92Nl,Nl,Z 

- Nl, ( - 9iNl,Zt^, + 91NIA) - V2\Nl,N;^,Zt^, 
+ 2V2KNl,Nl,Zi^, - Ar;3( - (^lAT^iZ.^i + ^sA^^^^tli 
+ v^A(Ar;4Z,-j3 + Nl.Zt^,))) (CIO) 

rfo^^o^^o^ = - ( - Z^i [Nt,^ [giNt.i - g2Nt,2) + 9iNt,iNt,3 - g2Nt,2Nt,z 

- V2X*Nt,,Nt,^ - V2X*Nt,4Nt,,) + 72^*^3 ( " 2«:*^ti5A^t25 
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A 
't-j,2 









XtiXt+^t3 

pi? 

Xt"iXtV*3 
XtiXt2^t3 






+ A* (iVi,3iVt24 + iV,,4iVt,3) ) + ^4 (A'tl4 (^liVt.l - 92Nt,2) 

+ ((7iiV,^i - (72iVt,2)iVM + V2\*{Nt,^Nt,^ + iV.^siVt.s))) (C.ll 



(C.12) 



(C.13) 

y^2V^t,lf/t.2^£ + V^*,2(^2f/t.l^Jl + A*f/,,2^2 



1 

71 



^2v;:if/;2^-2 + y:A92Ui,z^,, - xui^z^,^ 



■pij 



'^tlX^ji^ts 



pi? 






pL 

X?i 6*2^4*3 



^T/^^tSi [diNt^i - 92Nt22) 

3 

- 92Zt.;tlVt^l + 2_^ YeMji\3+ji^h2 



3 

/ .^lhhUR,t2.h^t3l 



(C.14) 
(C.15) 
(C.16) 
(C.17) 
(C.18) 

(C.19) 
(C.20) 



ii=i 






rre,* yiJ 

^L,t23lt:ih 



3 3 



't33+i2 



pi? 



Xti<^i2a2'^*3"3 



i2=iii=i 

3 



(c.2i; 



3 3 



ii=i 



J2 = lj\ = l 



(C.22) 






3 3 



i2=iii=i 



D 

t33+J2 



(C.23) 
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3 3 3 



133 O \ , . ^ . ^ 

(C.24) 
ii=i ji=i 

+ siv^ E E f^L;;,.>^«.i..^.^3+,.) (C.25) 

3 3 3 

r§^„,2-2«*3a3 " 3^°2."3 (2^C/1 2J ^£3+ilf^R,t2Jl^tll ~ 3 2^ 2_. ^«,ili2^£jl^i?,t2i2^*l^ 

ii=i i2=iji=i 

(C.26) 

rJ-.,.4 = -^^2f^;iEf^l:U^L. (C-27) 

ji=i 
rf. e,.;„ = ^ E E n:....^.^..f/^,...V^*.2 (C.28) 



Xtiet2'^t3 

i2=iii=i 



3 3 3 



^iV'2"2«*3<.3 ~*'^"2'°3 ( - 92UI1 2^ f^L',t2ii^2ii + f^t*2 E E ^L',t2il ^« .ili2 ^^3+^2 

ii=i i2=iii=i 

(C.29) 

r|.,.,.:3.3 = ^^«2,.3 E E ^4.2^.3..f^l*2.2^*.2 (C.30) 

i2=iii=i 



3 3 









"7| E ^^ftU2 E ^i,t2il^ejiJ2^t3l (C.31) 

i2=i ii=i 

, 3 3 

"7| E E ^^dij2 ^LMh ^kt232 ^t?i (C.32) 

J2 = lj\=l 

^ 3 3 

rgt^etjA?,^ ^ 77?J E ^^'.tli2 E ^i.'t2ii^e,JU2^t3l (C.33) 

' i2=i ji=i 

4 E E Knr2Uh.3J^h.nZti i^M) 

^^ J2 = ljl = l 

^~/^^ai,a2 2_^ ^R*tlj2 E ^L*t2Ji^d;hJ2Zt2i (C.35) 

^^ ^-2=1 ii = l 

-. 3 3 

-J=^ax,a2 2_^ 2-^ ^lhJ2^ LMh^ RMJ2^tH (C.36) 

^^ J2=lil=l 






^1 "^1 '-2*^2 '■S 



r^ = -i 
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-7^^aua2 2_^ UR,hJ2 Z^ f^L^aji^'^Jiia^tsl (C.37) 

V^ i2=l ii=l 

-. 3 3 

-pR — _x Sr^ Sr^ V* Tf^ Tf^ 7^ 

3 3 

^dt,o.iUt2C.2Hi^ " ^^*3i ^"l."2 2^ f^R',tlJ2 Z^ ^L,'t2il^'^Jli2 (C.39) 

72 = 1 ji = l 

3 3 
^Sicint^c.^^,- =^^it2*'^"l."2 Z^ Z^ ^«,iii2^L,tlii^i?,t2i2 (C-40) 

3 3 

pL — — 7__A \^ rr"'* \^ rr"'* V 7^ 

^ utia-i^ut^azhtg '■ /^""I'^z /_^ ^ RM32 Z_^ '^L,t2ii «JU2^t32 

i2=i ii=i 

1 3 3 

-7|'^ai,a2 2^ 2^ ^«,iii2^i,tiii^i?,t2j2^t32 (C.42) 

i2=iji=i 

, 3 3 

~/|'^Ql,a2 2^ ^ij',tU2 Z^ ^L,'t2Jl^«Jli2^t32 (C.43) 

i2=i ii=i 

-. 3 3 

- -7|5..«2 E E ^:...2f^L,.,.f/]^,2,2^4 (C.44) 

i2=iii=i 






^*lQl'"*2a2"^t3 



3 3 



^IQl "'i2"2 t3 

J2 = ij 



C.2 Two Fermion - one vector boson interaction 

'^X^v.^z^ = - \K,t2 [9130 + 92Ce) (C.45) 

T- w+ = - i^g2UP* , (C.46) 

r&,x?,7. = ^ (^i^e - 92se) {K^Nt^s - iV* 4iV*.4) (C.47) 

r|,x?,7. = - i (^i^e - ^2Se) (iV;3A^*23 - K,Nt,,) (C.48) 

r&,x?,^. = - i (^i^e + ^2Ce) (iV;3^^*x3 - N^.Nt,,^ (C.49) 

rfo^^c. ^^ = '- [gise + ^2ce) (iV;,3A^*23 - K^N^,,) (C.50) 

r|.^^^-H'+ = - l92{2V:^iNt,, + V2V,l,N,,,) (C.51) 

^4xr,< = - i^^ (2iV;2f/*2i - V2Nl,Ut,2) (C.52) 

^4x.-7. = ^(2^2V;:iSeV^*.i + V:^,^g^ce + ^2Se)v^*.2) (C.53) 

^&,x.-7. = ^(2^2f/;isef/t2i + f/;,2(^iC0 + 92Se)Ut,2) (C.54) 
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xt,xT,z, = 2 i '^92Vtli('eVt,i + V,l^ [ - g^se + g2Ce ) Vt,2 



^R 



xi^Xt^z^ 2 



25-2 f^tt 1^0^*21 + Ut^2 ( - 9iS0 + g2Ce ) Ut 



't22 






K,t2\9ice + 92Sej 






r^ =- 



-.i? 



^ti,t2(-fi'iS0 + fi'2Cej 
-^5'i5ti,t2Se 



:0ai,a2"ti,t2 



^licirftaaaTM ~ 3^lCe5ai,a2 5ti,t2 






5aua2K,t2 [^g2C0 + giSQ 






r 
r 



^ ofi'l^ai,a2^tl,t2'^© 



R 
L 



1 

'5ai,a2'^ti,t2 (35f2Se + fi'ice 



~ 6 

— :r 91(^0^01,02^11,12 



5ai,a2Sti,t2 [^g2Ce - gise) 



■^R 



2i 



rMtiaiMt2C2^M ~ r^9l^ai,a2Sti,t2Se 

C.3 Two Scalar - one vector boson interaction 



C.55) 
C.56) 

C.57) 
C.58) 
C.59) 
C.60) 
C.61) 
C.62) 
C.63) 
C.64) 

C.65) 

C.66) 
C.67) 
C.68) 
C.69) 



d-Haid* -y^ 






■p _ _ r 






i,a2 (( - ^92Se + gice) Yl ^£*^£i " 2^1^© Yl ^^S^n^Zz+n 

(C.70) 

3 3 

i,a2 ((3^2C© + (71 s©) Y Kn^Sn " 2(?iSe J^ Z,%^ ^£'3+,-^ ) 

(C.71) 
(C.72) 



ii=i 



• \ ^ ryU,* ryE 






ii=i 



ii=i 



(C.73) 
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1 ^ 

3 3 



(C.75) 



""tlQiftoQo^M 



(C.76) 

3 3 

3 3 

rg.^ej^z, = ^ (( - ^i^e + 92Ce) Yl ^St^Si - 2^iSe Y ^'^+3^+31^ 



ji=i ji=i 



(C.78) 

r.,^H,>- = - ^^2 (^.Ti^.ti - Z^.2ZU) (C.79) 

r^.,Ao^z, = ^ ( - ^i^e - 92Ce) {z^A - ^4^?2) (C.80) 

^A<i,^H+^W- = 292 [Zf^^Z^^^ + Z^^2Zt^2) (C.81) 

r^-H+z, =l{- 9ise + ^2ce) (^*ti*^*ti + ^*t2*^*t2) (C.83) 
C.4 One Scalar - two vector boson - interaction 

^h.^w^^w- = I9I {v,Z^,i + VuZ^^,) (C.84) 

^ht^z^z^ = 2 i9ise + 92Cef (vdZ"^ + VuZ^^^ (C.85) 

^H^^w+^^ = - 2^i^2 (t^d^itl* - VuZ^Jjce (C.86) 

^H,- w+z^ = \9192 (vdZtS - VuZt^^ye (C.87) 

C.5 Two Scalar - two vector boson interaction 



^ .2, 



^9t^w-v't^w+ -2^2K,t3 (C.88) 

^ 2 

rj>,^z,i>* z. = ^K,tz {9ise + 92Ce) (C.89) 
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i2=i 

3 3 

t33+J2 
J2 = l 32 = 1 

(c.9i; 



^ 3 3 



Th.,z.K.,z. = I {gise + ^2ce)' (^i^Ji + ^tTs^S) (C.93) 

^^0^1^-^°^^+ = 2^2 (^^t^l^tll + ^ti2^i32J (C.94) 

r^o^z.^o^z. = I (giSe + 92Cef {Zt^A + Zf^.Z^^,) (C.95) 

/fi;iy,-H+H^+ = 2^2 [Z^S^ti + ^4^2* ^4^2 j (<^-96) 

T/z-z./z+z. = ^ (-^ise + g2Cef {z^iZli + ^4^2* ^4^2) (C-97) 



r^- 



C.6 Four Scalar 

We define 



3 33 333 

^2 = - ^( XI ^S+i3 Yl Yl ^f,hh^f'nJ2Z[,3+J2 + YY1 ^S Yl ^Ihh^f^hhZ!^ 



333 33 

F 
33 

33 = 1 J2 = lji=l J3 = 1J2 = 1 il = l 

(C.99) 
3 3 



^3 =1 (a E <3%,. E ^/:...2<.. + A* E E <^* ^/....<3+,.) (c.ioo) 

i2=i ji=i i2=iii=i 

^4 = - -^ ( E <3V.2 E ^;.U2<. + E E <;;^/....<3...) (cioi) 

J2 = l il = l J2 = ljl = l 

^ =1^ ((^A^i + ^Cl9l) E <^* <.. + 2^? E <'^..<3..0 (C.102) 

^6 =^ (a e <3%,2 E ^/:...2<.. + A* E E <^* ^/.i..<3+,.) (C.103) 

i2=i ii=i i2=iii=i 

3 33 

^7 = - ^( E ^tl3+i3 E E ^/Jli3^/Jli2 ^t^3+i2 
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3 3 



i3=ii2=i ii=i 

3 3 3 



t2 J3 



i2=i ii=i J2=iii=i 

-.3 3 3 3 

^9 =2 ( - ^ ZZ ^S+i2 ZZ ^Ihh^lh ^^^ZIZI ^5i^/,iii2^, 



t2 3+J2 



t2 3+i2 



(C.104) 
(C.105) 
(C.106) 



i2=i ji=i 



With this definitions often appearing terms in the vertices involving squarks and 
sleptons are given by 

A =A, with Yf ^ Yd,Tf -^ Td,Z^ -^ Z'',Cl ^ l,C]j ^ l,Ci ^ 1 (C.107) 

Ui =A, with Yf -^ Yu, Tf -^ T„, Z^ ^ Z^, Ci ^ 1, C]j ^ -2, Cl -^ -1 (C.108) 
Ei =Ai with Fj -^ Fe, Tf -^ Te, Z^ ^ Z^, Cl ^ 3, C]j ^ -3, C| -^ -1 (C.109) 



r- 


0q,-^^q,2 




+ /^3 


^dt,^^A^,^dl,^^A^l^ 


'^ai,a2 




+ D^ 


r«t-^ci/lt2«t3«3H 


^ 0q,-^^q,2 




+ f/3 



+ f/3 



-Dl(^^t2l^t4l ~ ^t22^t42J + -02^t2l^t4l 



^H vi? 



7// ryH 



^t22^t43 + ^i23^t42 



-^l(^^t2l^t4l ~ ^t22^t42J + -^2^*21^*41 



rvA rvA _ rvA ry A 
^t22^t43 ^t23^t42 



H ryH 

ti2 



Ui^y^tii^ui ~ ^t22^i 



ryH ryH . ryH ryH 
^t2l^t43 + ^t23^t4l 



^2^t22^t42 



(C.llO) 



(C.lll) 



(C.112) 



rA ryA 
'U2 



^2 ^i2 2^*4 2 



(C.113) 

H ryH 

t42 



^l(^t2l^t4l ~ ^t22^) 

_ ryA ryA _ ryA ryA 
^t2l^t43 ^t23^t4l 

^ et^ht^el^ht^ =Eiy- Z^^^Z^^-y + ^4^2^*42] + -^2^*21^*41 + -^3 (^^t22^t43 + ^t23^i 

(C.114) 

T _ -p ( ryA ryA , ryA ryA \ \ T^ ryA ryA , Tp ( ryA ryA ryA ry 

et^Ah,t2e;^Ah,t4 -J^iy- ^t2l^t4l + ^t22^t42J + -^2^t2l^t4l + -^^3^^ - ^t22^t43 " ^t23^i 

(C.115) 

^ht^ht^ht^ht^ = 7(^*11(^*21 ( ~ 35f+^t3l^t4l + '^^t32^t42 ~ MM ^t33^t43 



A ryA 

t42 



+ Z^^, [\ZZ2ZZ1 + A^t^i^.^2 + 4Re{«:A}Zj3Z,^; 



+ 2-^*23 ( ^'^* ( ^t32^t43 + ^(33^*42 



+ A* ( [kZ^^, - 2AZji) ZJ3 + ZJ3 [kZZ^ - 2AZji) ) ) ) 
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7H 

^t42 



+ 2^423 (^A/t* (^Zj^iZj^3 + Z^,^^Z^^]^j 

+ A* ( (/.Z,^, - 2AZJ2) ^,^3 + ^£ {>^ZZ, - 2XZ^, 

+ 2Z/f3 (a* [zZ^ [Zl, [kZZ^ - 2\ZZ,) + ^^2 ('^^t?! - 2AZ, 

+ Z,^i ( («:Zj2 - 2AZJ,) ZJ3 + ZJ3 [kZ^^, - 2AZji) ) 

+ Z^^, ( («:ZJ, - 2AZJ2) ZZ, + Z,^3 (/.Z,^, - 2\Z^^,) ) ) 

+ /.* (Z,^3 ( - 12/.Zj3Z,^3 + \Z^^,Z^^, + AZ,^2^,^i) 

+ -^(^^ta 1(^^*32^*43 + ^433^*42 j + Zt22\Z^^^Z^^^ + Z^^^Z^^^ 

(C.116) 

T/i.^htjAOgAO^ = ^(^432(^*42 (A^til^tal ~ 4|ApZ^^3Z^23 " fi'+^ti2^t22J 

+ 2(«:A* + \K*) ( - Z,^iZ,f3Z,^3 + Zt,,[z^,A + ^tfs^t^i 

+ ^tgl y " ^t4l (^ ~ A^ti2^t22 + 4|A| Z^^3Z^23 + ^f^Z^^^Z^^l 

+ 2 (/.A* + A/.*) ( - Z^^^ZZA + ^*13 (^?2^& + ^?3^, 

+ 2^4 (a* ( - Zf^, {Z^^, {2\Z^^, + kZ^^,) + Z,^2 (2A^,^2 + ^^^t^i) ) 

+ '^[Zt^l (^^ti2^t23 + ^ti3^t22J + ^t42 (^^til^t23 + ^ti3^t2l 

+ /.*(- Z,-^3 {AkZ^^.Z^^, + AZ,fiZ,f2 + A^?2^&) 

+ A( ^i4l( ^ii2^t23 + ^ti3^t22 ) + ^t42( ^til^t23 + ^ti3^t2l 

(C.117) 

+ Z,^2 (A^£^i + ^+<2^i) ) - ^S* (4AZ,^3^t?3 {^*Zl, + A*Z+i 

+ '^ti2(^AZ^3iZj^2 + 5'-^t32^t4lj + ^til(^AZ^g2^t42 + fi'+^t3l^t4lj j 

(C.118 

+ Z,^2 (a^4^m + A^i^i + 4Re{«:A}z4z,^3) 

+ 2Z4 (A/.* (^4^4 + Z^l.Zf^,) + A* ( ( - 2AZ,^i + KZfl,) Z, 



't22 
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+ Zf^,[-2\Zf^, + KZ, 



A 

t4,2 



+ ^ti2 ( ^t22 






35'+^t32^t42 + "^^tzl^til - 4|A| Zjg3Zj^3 

+ 2^t23 ('^'** (^t3l^t43 + ^t33^t4l j 

+ A* ( ( - 2AZ4 + ^Zt,i) Zt, + ^4 ('^^tli - 2AZ, 
+ 2Z,^3 (a* (z4 (z,^i ( - 2\Zf^, + /.Z,^^) + ^4 {^Zf^i - 2A^, 
+ Z,^i ( ( - 2\Z^^, + /.Z^) Z,-j3 + Z,^3 ( - 2\Zt^, + «:Z,^2 
+ Z4 ( ( - 2AZ4 + ^Z^^i) Zt, + ^,^3 ( - 2A^,l2 + '^^ili 
+ K* (z4 ( - 12/.Z4Z,^3 + AZ^iZ^ + A^4^i) 

+ A(^i2 1(^*32^*43 + ^t33^t42 ) + ^t22( ^t3l^t43 + ^t33^t4l 



A 

t4,2 



i* ryA 



-.2 ryA 



(C.119) 
^?i^t2^t3-f^t!i ^ 4 V" ^*2^* (^4A*Zj^3Z;^3 (^AZ^^2 - /^^tti j + ^tti (^fi'-^tti^tt2 ~ AZjg2^t4 
+ Z,^2(^+^4^i - A^ili^i)) - ^Sf (4AZ,-^3^43(^*^i - '^*^, 

^(42 







* '7 + 
(42 



+ ^ti2(5'-^t32^t4l ~ ^^t3l^t42 ) + ^til(fi'+^t3l^t4l ~ ^^t32^t42 



(C.120) 



+ ^til ( ~ '^9+Zt2l Zp^lZt^l + A^(22 ( ^t3l^t42 + ^t32^t4l 



dtiaii>t2dl^a-i'^t^ 24 



(C.121) 



,1 
'5ai,a3 (5*2, t4 (2^? 5^ Z5*+j.^^g3+ii + (3^2 + 



ii=i 



ii=i 



+ 5,,,, {2gl J2 Z,%^Xs.n + {^9l + ol) E <i:^4)) 



i2=i 



J2 = l 



dtiaiet2dt^a-^'^t4 



r)^"ai."3 



2^72 V z^'* 



/ . ^t23+il^t43+il \^^ /^ ^tl3+J2^t33+J5 

ii=i 



i2=i 



32 



3 

i2=i 



(C.122) 

rvD,* yD \ 



6 6 

^t2h^Ujl (2fi'l 2^ ^tl3*+J2^23+J2 + [9l - 35-2 j 2^ ^hjl^t!ij2 



jl=l 



J2 = l 



J2 = l 



29 



2g] Y. ZX.^Sw. + («? - 39O Z ZffizS.) E Z£2.' 



3 

ii=i ii=i i2=i 

3 33 

_L 0^2/0 V^ 7D,* 7D , \^ ryD,* rvD \ ST^ yE,* rvE 

+ ^9\ y- Z^ ^tl3+ii^t33+il + Z^ ^tlil^tail j Z^ ^t23+J2^t43+i2 

ii=i ii=i i2=i 

33 33 



+ 24(^ 2^ 2^ Z^2ii^e,iiJ2^t43+i2 Z^ ^ti3+J4 Z^ ^IjsJA^hj: 

i2=iji=i i4=i i3=i 

3 3 3 3 

+ ^ ^ ^tlii^'i.iii2^t33+i2 Z^ ^t23+i4 Z^ ^eJijiZt^j^Jj 



J2 = lji = l J4 = l J3 = l 

(C.123) 
3 3 

3 33 

- 12 Z^ ^tl3+i3 Z^ Z^ ^djii3^rfjli2^t33+J2J^t4l 

+ 12 E E ^5: E ^:.3..^«.2..^4) ^i) (C.124) 

i3=ii2=i ii=i 

r,,^,,^,^,;^ = - 1 (^(^2 ^ (^2) (2<5ti,t3^t2,t4 + 2St,,tA.,t,) (C.125) 

3 3 

^i>t,ut2C2'>t,ut^o.4 = ^^a2,«4 ('^ti,t3 (( - 35-2 + alj E ^42^*1 ^i«l ~ ^df E ^t23+ii^tl3+ii 

+ ^..3 (( - 3..^ + 9l) E <^2^.^.2 - ^9! E <3%.2^1l3..2)) 

i2 = l J2 = l 

(C.126) 
3 3 

^ut,et2%et^ = - ■^[K,t3[[9l - 9lj E ^t2h^Uh ~ 2fi'l E ^t23+ii^t^3+iij 

33 33 

I „2 / \ ^ ^S,* rvu \ ^ rvU,* rvE , \ ^ ryV,* ryE \ ^ ^S,* ^1^ A 

+ 5-2 1^ Z^ ^Hh^Hh Z^ ^tli2^t4i2 + Z^ ^tlil^t4il Z^ ^Hh^Hn) 

ii=i i2=i ji=i J2=i 

3 3 3 3 

+ ^ E E ^tlil^e'ili2^t43+J2 E ^t23+i4 E ■^e.i3i4^t3i3 j 
J2 = ljl = l i4 = l J3 = l 

(C.127) 
^ vt^ht^vtht^ = ~ 4 (fi'i + fi'2 j^ti,t3 y^tii^ui - Zti2^u2) (C.128) 
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^ut.AO^ul^AO^ = ~ 4 (fi'l + 9l)Sti,t3 (^t1l^t!l - Zf22^U2) (C.129) 

3 3 3 

j:i=lJ2=l ii = l 

+ [9l-9l)z^j6,,,,Zl,) (C.130) 

3 3 3 

r - J-A r _ 4n2 \^ 7*^'* 7*^ r \^ 7-^'* 7-E _ 9 "SP yE,* yE 

ii=i i2=i i2=i 

Jl=l J2=l i2 = l 

+ ((39? + 9?) E ^S^ln - 29? E ^5+.^.!3+,.) E <J^W. 

ii=i ii=i i2=i 

3 3 3 

ii=i ii=i i2=i 

(C.131) 

3 3 

3 3 3 

+ 3 2^ 2^ ^tijt. 2^ yrf*j3j-^Fdj2jiZi3j3J jz+i 

i3=ii2=i ii=i 

- ZS{{3d + 9?) E ^a^g., - "9? E ^X.^fi+i, 

+ 12 E <3%,3 E E Knny^^nnZl!,s^n)K2) (C.132) 

i3=i i2=iji=i 

3 3 3 3 

r — —-( n'^ \^ 7^'* 7^ \^ 7-^'* 7-^ _L /,2 Y^ 7^,* ^B V^ ryE,* r^E 

J- 6^16*26*36*^ - g l^i/l Z^ ^tlil^i4il Z^ ^t2i2^i3j2 "^ ^2 Z^ ^tin^Un Z^ ^t2J2^t3J2 

ii=i i2=i ii=i i2=i 

3 3 

_ n 2 ST^ yE,* yE ST^ yE,* yE 

^yi Z^ ^tl3+ji^t43+il Z^ ^t2i2^t3i2 

ii=i i2=i 

3 3 3 

_ 9/,2 \^ 7^'* vE ( \^ yE,* yE _ 9 V^ yE,* yE 

^i/l Z^ ^t23+jl^t43+il \^ Z^ ^tiJ2^t3J2 ^ Z^ ^tl3+J2^t33+i2 

ii=i i2=i i2=i 

+ E ^fi^£. ((»? + ^0 E <S^4 - 29? E ^5«.^.< 
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0„2/ V^ ryE,* ryE ST^ 7^,* rjE , V^ n'E,*ryE ST^ ^E,* rjE 



■i 6 



ii=i i2=i ji=i i2=i 

3 3 

't33+J2 

3 3 

I \^ yE,* ryE V^ rvE,* ryE 

+ Z^ ^t23+il^t33+il Z^ ^tl3+J2^t43+J2 
il = l 72 = 1 

3 3 3 3 

, „2 V^ ^S,* rvE Sr^ rvE,* rvE , 2 ST^ yE,* yE \~^ yE,* yE 

+ 9l 2_^ ^t2h^t3n Z^ ^hJ2^UJ2 +92 2^ ^t2h^hn Z^ ^hJ2^UJ2 
jl=l J2 = l jl = l J2 = l 

3 3 3 3 

0^2 V^ yE,* yE \^ yE,* yE , „2 V^ yE,* yE \^ yE,* yE 

- ^9l Z^ ^t23+jl^t33+il Z^ ^tli2^t4i2 + fi'l Z^ ^tlil^tail Z^ ^t2i2^t4i2 

ii=i i2=i ii=i J2=i 

3 3 3 3 

+ fi'2 Z^ ^tljl^tail Z^ ^t2i2^t4i2 " ^5'l Z^ ^tl3+ii^t33+jl Z^ ^t2i2^t4i2 
il=l 32 = 1 jri=l J2 = l 

3 3 3 

9 „2 / V^ 7-B,* r7_B _ r) \^ yE,* yE \ \^ yE,* yE 

^9l y 2-^ ^tlji ^t'ih ^ Z^ ^tl3+il ^t33+jl y Z^ ^t23+J2 ^t43+J2 

ii=i ii=i i2=i 

3 3 3 3 

7E 

32 = 1 j\ = l J4 = l J3 = l 

3 3 3 3 

, V^V^^B,*^ 7^ \^ yE,* X^v-* 7-^ 

^ 2^ 2^ ^tin^e,jiJ2^t43+J2 Z^ ^t23+i4 Z^ -^e,i3J4^t3i3 

i2=iii=i j4=i i3=i 

3 3 3 3 

I Y^ V^ yE,*^ yE ST 7^'* Y^ V* 7^ 

^ 2^ 2-^ ^t2Jl^e,jiJ2^t33+J2 2^ ^tl3+i4 2^ ^e,J3J4^t4J3 
J2 = ljl=l 74 = 1 J3 = l 

3 3 3 3 

I \^\^7S,*V 7'E' \^ yE,* \^ v* yE 

"^ Z^ Z^ ^tlil e,jii2^t33+i2 Z^ ^t23+J4 Z^ -^e,i3i4^t4J3 
i2 = lil = l i4 = l i3 = l 

(C.133) 

3 3 

^t2i[[9i + 92 J Z^ ^hh^tih ~ '^9f 2^ ^tS+h^t^iZ+h 



6 6 6 6 

^'^\Z_^Z_^ ^t2Ji^ejiJ2^t43+J2 Z^ ^ti3+i4 Z_^ ^e,J3J4^t: 






3 3 



+ ^ Z^ ^tl3+i3 Z^ Z^ ■^ejijr3-^e,jriJ2^i33+J2J^i4l 

i3=i i2=iii=i 

+ K'2 ( - 2.? E <c...^.f3... + {9I + .2^ E ^S^^O^*^^ 

ii=i ii=i 

(C.134) 
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C.7 Three Scalar 



H 



(C.135) 

(C.137) 
r^...i'i.*,.,A03 =5,„.,(f/8^4 + U,(v,Z^^, + VsZ^^,)) (C.138) 

(C.139) 

Te-.^e-^AO^ = i?8^,i + E, (vsZf^, + V.Z^^,) (C.140) 

3 



+ V2 (.„ E <3%,3 E E Kn^J^^nnZ^ 

i3 = l 32 = 1 jl=l 

J3 = lJ2 = l il=l jl = l 

- 2(2 E <3V,. E ^«v..^£. + ^(-^^* E E <;;^^....^*' 



i3=ii2=i ii=i ji=i 

33 33 

7D 

't23+J2 

i2=i ii=i i2=iii=i 

3 33 

+ Vrf 2^ ^il3+J3 E E ^«*ili3^c'jli2^t23+i2 
i3 = l J2 = ljl = l 

+ -« E E <i2 E ^:m.>^«....^4)))) (c-141) 

i3=ii2=i ji=i 

r..,.4..3 = -^(^? + ^2)'^*i,*2(^<^^S-^.^£) (C.142) 

r.,.*,Hr3 = i (^^i* ( - ^'"" ^ ^*-^"- + 2..A* E E Z^;^J.nnZ. 

ji=i i2=iii=i 

+ zts ( - 72,,^, E Zt;*Xn + 4 E E ^r4Te,,,.^.f3-... 
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ji=i 



33=132=1 



E 

t2J3 



(C.143) 



+ Z^.z [V22Re{T^]zZ, + 2A/.* (t;,^^^ + v^Z^^,) 

+ 2A* (A3ZJ3 - 2vs\Z^^, + t;s«:^£))) 

+ Z^^ (^£ ( - 3^^„^,f2 - 4t;.| Al^Zjg + t;,AZ,fi 

+ Z^^, (t;„AZji + t;,AZ,^2 + Ai^.^a) 

+ Z,^3(v^Re{T4Z,^, + 2A/.*(t;,Zj3 + v.Z^^, 

+ 2A* ( - 2t;,AZ,f2 + A2^t?3 + ^s«:^Ji) ) ) 

+ Z,f3 (v^( - 4Re{T4Z,f3Z,f3 + Tl {z^^.Z^, 

+ T, [Z^,,ZZ, + ^,^2^£) ) + 2k* (aZ,^2 (^.^S + ^.^£ 
+ A^.^i (t;.^£ + ^n^&) + ^t^3 ( - ISt^s/t^Ja + Vd\Z^,2 + ^«AZji 
+ 2A* (Zj3 (AsZ,^! + A2ZJ2) + ^,?i (Aa^Ja " 2t;.AZji + VskZ^ 
+ Z,f2( - 2t;,AZ,f2 + A2ZJ3 + t;.«:Z,fi)))) (C.144) 

. Z,^, ( - 4t;,Re{A/.}Z,^3Z,f2 + v^2 RelT^Z^Z,^^ 
+ 4t;,Re{A«:}z4z,^3 + v^2Re{T,}z4z,^3 - AvJie{\K]zf^,Z^^, 






^^H^lAz ~ I 



H 



vH 



+ Z^AAv,\\\'Zii, + givaZ^^, - v^XZ^^ 



+ Zf^,[2v,K\*Zt,,Z^^, - V2T:Z^^,ZI;, - V2T^Zt,,Z^^, - 2vsk\* Z^^.Z^,, 

- V2TlZt^,Z^^, - V2nzt^,Z^^, + 2v,K\*Zt^,Z^^, 

+ Z4 ( - 4t;.| ApZ/f3 - glu^Z^^, + t;,AZ,fi) 

+ 2\K* ( - t;.Z,-JiZ,f3 + Z,^3 (t;,Z,f3 + t;,Z,^,) ) ) 

+ Z4 ( - V2 ( - 2Re{T4z4z,f3 + T,* (z.^^Z.f^ + Z^Z.fi) 

+ T, (z,^iZ/f2 + ^4<i) ) + 2A* ( - Z4 ( (st^^A + Vuk) Z^^, 

+ 2k* (az4 (^rf^/fs + VsZ^^,) + AZ,^, (t;,Z,^2 + ^«^5 



"^n^tiS 
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- Zf^,(^AvsKZ^^, + v,\Zl^ + t;„AZ,fi)))) (C.145) 
r.,^^+^- =\{- ZtS [K^ [f-VuZl, + v,\Zl^ + Zl, {glv.Zl, + vJ.Zl,) 

- Zi; (Z,f, (^^d^,t2 + ^«AZ+i) + ^,^2 (^^«^tt2 + vAZl,) 

+ 2Z,f3(2t;,A*(«:Z+i + AZ+2) +v^T;Z+i))) (C.146) 

- Zjf (2 {2v.M* - V2T,) Zf^, + vAZf^, + v^XZf^,) Z+2) 

(C.147) 

D. One-loop tadpoles 

In this and the subsequent Apps., particles that are denoted with a hat, e.g. /ij, are 
the unrotated external states. In the corresponding vertices the associated mixing 
matrix has to be replaced by the identity matrix. Moreover, we have summed her 
and in the subsequent section in all the vertices implicitly over the colour indices of 
quarks and squarks. 

At the one-loop level, the expressions for the tadpoles of eq. (2.13) are given by 
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E. One-loop self-energies 

The definitions of tlie scalar one-loop functions and their explicit analytic expressions 
can be found in ref. [14]. 

E.l Self energy of Z-boson 

In agreement with ref. [13] we obtain for the transverse self-energy of the Z-boson 
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E.8 Self-energy of sneutrinos 
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